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Recently, several simple quantum mechanical toy models of black hole evaporation have 
appeared in the literature attempting to illuminate the black hole information paradox. We 
present a general class of models that is large enough to describe both unitary and nonuni- 
tary evaporation, and study a few specific examples to clarify some potential confusions 
regarding recent results. We also generalize Mathur's bound on small corrections to black 
hole dynamics. Conclusions are then drawn about the requirements for unitary evapora- 
tion of black holes in this class of models. We present a one-parameter family of models 
that continuously deforms nonunitary Hawking evaporation into a unitary process. The 
required deformation is large. 



I. INTRODUCTION 

Hawking's calculation of black hole evaporation [lj leads to a direct conflict between general 
relativity and the unitary evolution of quantum mechanics (2). This prompted the suggestion that 
the requirement of unitary evolution should be relaxed [3, 4], and pure states allowed to evolve 
into mixed states. Such nonunitary evolution, however, seems problematic O [j] 

Moreover, since black hole evaporation is a very slow process involving a large number of 
emitted particles, and since one only expects to start recovering information after about one half 
of the radiation has been emitted [8], one might imagine that unitarity is restored by the accumu- 
lation of small (say, nonperturbative) corrections over the course of the entire evolution process, 
cf. 13. To evaluate this claim, Mathur [10] introduces a qubit model of black hole evaporation 
and then derives bounds on the entanglement entropy of the emitted radiation, which show that 
this scenario is not possible. Below we generalize his result to consider more general deformations 
of the pair creation dynamics. In particular, for small changes in the evaporation process the en- 
tanglement entropy continues to increase, therefore (barring remnants and variants thereof) the 
evolution is not unitary. 

Thus, for unitarity to be restored, one needs to make large corrections to the semi-classical 
evaporation process described by Hawking. Recently, several specific unitary models have been 
introduced fnl - [i"3| as proposed alternatives to Hawking's semiclassical evolution. As discussed 
in [14|, models of this kind could be called "burning paper" models that involve large corrections 
to the Hawking evolution. While these models are unitary, they are typically written in a way 
that makes it difficult to compare to the semiclassical evolution and to see how Mathur's bound 
operates. The difficulty arises because the bound is derived in terms of dynamics in an ever- 
enlarging Hilbert space, whereas unitary models are typically written as dynamics in a fixed- 
dimensional Hilbert space. 

The primary goal of this paper is to clarify the meaning of Mathur's bound in pO} . In partic- 
ular, until one writes unitary models as (large) corrections to nonunitary, semiclassical evolution, 
the meaning of Mathur's result remains obscure. Previous investigations of corrections in this 
context either directly introduce unitary models that are difficult to compare to Hawking evap- 
oration ITTI - IT3| ; or consider small corrections to Hawking evaporation that do not produce unitary 
evolution when made sufficiently large |T3j [25f, and thus are unconvincing illustrations of the result 
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in [TUJ. A second goal is to characterize what kinds of corrections produce the desired unitary evo- 
lution. That the corrections need to be large is a necessary condition derived in llTOl , but sufficient 
conditions have not been discussed in the same way. 

In order to address the above points, it is necessary to introduce a general model space that 
provides a uniform language to discuss both unitary and nonunitary black hole evaporation. This 
allows us, for example, to continuously deform the semiclassical Hawking evolution to unitary 
evolution. One can then explicitly see that the deformation is large in an appropriate sense, and 
therefore in agreement with a suitable generalization of Mathur's argument. Let us emphasize 
that it is not our intention here to advocate for nonunitary evolution, only to demonstrate that 
unitarity demands there be a significant alteration of the traditional semiclassical evolution. To 
restrict ourselves to unitary evolution at this stage would be to beg the question. 

In Section |nj we introduce a very general framework for qubit models of black hole evapora- 
tion that is appropriate for both unitary and nonunitary evolution. Most of the discussion focuses 
on how to interpret the models, and their connection to ideas in quantum information theory. 



In Section III we apply the formalism to a number of sample models; some of the models were 
chosen because they were discussed previously in the literature, and others because they illus- 
trate some interesting issues. In Section |IV) we briefly review Mathur's argument against small 
corrections restoring unitary evolution, and generalize the main entanglement entropy bound to 
allow arbitary deformations. Mathur's original result IfTUl only explicitly considers one kind of 
perturbation. In Section|Vj using observations from Section III we discuss what conditions ensure 



unitary evolution. In Section [VT] we present a one-parameter family of models that continuously 
connects Hawking evaporation to a unitary model; one sees clearly that the deformation required 



is large. In Section VII we conclude with some brief comments. 



II. THE GENERAL MODEL 

Before presenting our model, we make a few preliminary comments on nonunitary evolution 



in Section II. 1 Then in Section II.2 using some results from quantum information theory, we 
explain how to describe nonunitary evolution that still has a good probabilistic interpretation. 
Along the way, we clarify some potential confusions regarding previous work. Finally, we present 



our general class of models in Section II.3[ discussing the physical interpretation in Section II.4 



ILL Nonunitarity 

When we model the evolution of a closed quantum system, the state of the system is given by 
a ket \ip(t)) that satisfies 

\m) = u(t)\m) (2-i) 

for a unitary time evolution operator U (t). We may equivalently write the state of the system as a 
density matrix p(t) = \ip(t)) (ip(t) \ where p(t) satisfies 

p(t) = Up(0)Ul (2.2) 

Unitary evolution satisfies several nice conditions, namely^] 
1. Linearity 



2 We do not claim that these are completely independent. 
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2. Preservation of the norm: unit norm states evolve to unit norm states, ensuring that a 
probabilistic interpretation makes sense. For density matrices, the desired condition is that 
unit-trace, completely positive density matrices evolve to unit-trace, completely positive 
density matrices. 

3. Invertibility: previous states can be found from the current state. 

4. Purity: pure states evolve to pure states. 

However, it has sometimes been suggested (2]-31 that theories of quantum gravity (especially 
in the presence of black holes) will not be unitary. Let us note that the negation of unitary is 
ambiguous, since it is not clear which of the above conditions is relaxed. Let us consider three 
illustrative mappings. 

1. This evolution does not conserve probability — the norm is not conserved; however, pure 
states still evolve to pure states, and the evolution is invertible: 

m »\\i>). (2.3) 

This kind of evolution is sometimes useful in modelling a system that decays into some- 
thing that is outside the model. For instance, in modeling alpha decay. An equivalent 
way to describe this evolution is to say that the system has energies with an imaginary 
part. For a fundamental description that includes all degrees of freedom, however, the 
nonconservation of probability is nonsensical. 

2. This evolution preserves the norm, is invertible, but evolves pure states to mixed states: 

IV'i) ^ Pi = 2 1^0 (^il+j l^ 1 ) l^ 2 ) ^ p2 = \ 1^ + \ 1^2) (<h\ ■ ■ ■ (2-4) 

While this model evolves pure states to mixed states, information is still preserved (in a 
weak sense)J^]and probability is conserved. 

3. This evolution conserves probability, evolves pure states to pure states, but is not invert- 
ible. 

In this model one cannot reconstruct the past from the current state, and therefore infor- 



mation is not preserved. In Section III.4 we introduce some models of this type. 

In the black hole information paradox, one considers some initial configuration of matter \ip m ) 
that collapses into a black hole, and then completely evaporates to radiation in a thermal mixed 
state p ~ e~P H . Since the final state is both mixed and independent of the initial state, the evolu- 
tion both fails to be invertibile and pure in the above senses. 

Suppose for the nonce that Hawking's original argument is correct [2], and the fundamental 
theory of quantum gravity is not unitary. We can no longer write evolution as in ( |2.2| . We restrict 
our considerations to evolution that satisfies conditions [I] and |2j but not necessarily condition [3] 
and|4j Then, assuming some very basic conditions that ensure a good probabilistic interpretation, 
we can write the most general possible evolution in the operator-sum representation (cf. |16|) 

p(t) = 2>*p(0)4 (2.6) 



3 Provided large ensembles of identical copies of the system, one can with some confidence distinguish distinct density 
matrices; however, this usage of the phrase "information preservation" is not canonical, and is problematic if one 
starts considering density matrices which are very close to each other. This would mean in an experiment repeated 
many times with identical initial conditions, one could reconstruct some information about the initial state from 
the final state. In a technical sense, however, quantum information is lost. We make the distinction here, since the 
semiclassical description implies that information is not preserved even in this weak sense. 



4 



for some set of operator Ef. that satisfy the completeness relation 

Y,E k {tfE k (t) = I. (2.7) 

k 

This is one way to write the evolution of an open quantum system, and the transformation from 
p(0) to /?(/) is called a "quantum operation" in the quantum information theory literature [16|. 
The operators E k determine the evolution of the density matrix. When there is only one E k , the 
evolution is unitarvFl 



II.2. Quantum Operations 



As in |T0l[T2U l5l, we model the Hawking evaporation process as a discrete set of mappings on 
qubits. In the initial state, the system consists entirely of matter in a pure state. There have been 
some suggestions [17] in this context that the entanglement between the initial black hole-forming 
matter and the outside matter plays an important role; we do not address these issues at this time. 
The initial state is modeled as a set of n "matter qubits": 

Po = l^o) {ipo\ l^o) G span{|gig 2 • • ■ q n )} , (2.8) 

where each q is a qubit, a quantum state labeled by or 1. After a sequence of intermediate steps, 
the end state consists entirely of radiation (again, we are assuming no remnants), modeled as a 
(possibly mixed) density matrix acting on n "radiation" qubits, pf. Throughout the evolution, we 
keep the total dimension of the Hilbert space fixed. This is certainly true for the unitary evolution 
of closed systems, but here we put it in as a reasonable assumption. We are motivated in part by 
the black hole's entropy. The black hole initially has S ~ M 2 , which entirely radiates away on 
the time scale ~ M 3 with an emission every ~ M. This is consistent with a model having a fixed 
number of physical qubits. 

Following [T2| , we use hats to distinguish the internal black hole qubits from the external 
radiation qubits. The hatted qubits represent all degrees of freedom that are inaccessible outside 
the black hole; unlike flOj, we do not distinguish between degrees of freedom from the initial 
matter, from gravitational interactions, or any that arise during the evaporation process. We write 
basis elements for the final state as 



{\q n qn-i---qi}}, (2.9) 

where we have put the labels on the qubits in reverse order for reasons which should become 
clear. At the iih step, we have a density matrix acting onn-i black hole qubits and i radiation 
qubits, so that the total dimension of the Hilbert space is fixed. The evaporation concludes on the 
nth step when there is only radiation: 

Po -> Pi -> > Pn-i -> Pn (2.10) 

In general the mapping from po to p n = pf should be a quantum operation, which is the com- 
position of n quantum operations (one for each emission). Therefore the total evolution from po to 



p n (and each intermediate step) may be written in terms of some E^s, like in Equation ( |2.6[ ). Quan- 



tum operations, however, may be written in an equivalent, alternative form that is more natural 



4 Note that we model evolution with discrete time evolution, and do not discuss continuous evolution, which might 
be governed by the Lindblad equation. 
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when discussing black hole evaporation. This form connects more directly with the discussion 
in IflOl . 

We motivate this alternate form, by noting that any mixed density matrix may be "purified" 
by enlarging the Hilbert space. For example, a density matrix of the form 

p = Pl \A)(A\+ P2 \B)(B\ (2.11) 

with orthonormal {\A) , \B)} can be purified by introducing the orthonormal states {\a) , |/3)}, and 
defining 

|*) = \A) <g> \a) + Vm\B)<E) |/3) . (2.12) 

Then, one sees that p is the reduced density matrix found by tracing out the new degrees of free- 
dom. Let us emphasize purification is a formal mathematical operation, and not a dynamical process. In 
particular, for us, the new kets that we direct producted into the Hilbert space do not correspond 
to any physical degrees of freedom^] Roughly speaking, then, we can imagine purifying each of 
the piS by enlarging the Hilbert space, and then the evolution in this enlarged Hilbert space would 
be unitary. 

As it turns out, any quantum operation from say po to p n may be written in the following way: 

Pn = tr aux [C/(paux ® Po)U j ], (2.13) 

for a unitary transformation U acting on some auxiliary degrees of freedom as well as the phys- 
ical degrees of freedom. In particular, if po acts on a d-dimensional Hilbert space, then we need 
introduce at most a ci 2 -dimensional auxiliary Hilbert space to write the most general quantum 
operation in this form |[l6fl. Thus for our n-qubit system, we need only introduce 2n auxiliary 
qubits to capture the most general evolution of density matrices. 
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(c) final 



FIG. 1. Here we illustrate the black hole evolution with auxiliary degrees of freedom. At first we have 
only the initial black hole degrees of freedom represented by dots. At later stages in the evolution, we 
introduce some new degrees of freedom that corresponding to the infalling negative energy particles in the 
Hawking process (squares), as well as outgoing radiation (wavy lines). We have thus increased the size 
of the Hilbert space, but this should be thought of only as a convenient way to parametrize potentially 
nonunitary evolution. To return to a fixed dimension model, one needs to trace out auxiliary degrees of 
freedom. In general, at intermediate steps, it is not clear which degrees of freedom should be thought of as 
auxiliary, so we illustrate just one possibility. On the other hand, it is unambiguous that in the final state all 
of the black hole degrees of freedom (dots and squares) are auxiliary if there are no remnants. 



In this language, we can think about the semiclassical Hawking evolution in the following 
way. We start with n black hole qubits initially in a pure state. We imagine that n might be 



5 If one does treat the new degrees of freedom as physical in the final state, then one is considering a remnant scenario. 
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roughly given by the entropy of the black hole|j At the first time step, a pair of qubits is created 
at the horizon in an entangled state, 

^(|6>|0> + |i)|l». (2.14) 

The zero represents no particle and the one represents a particle. We refer the reader to [lOl [15| 
for a thorough discussion on the origin of this description; see also |(T2|. We have now added 
two new qubits to the system, increasing the size of our Hilbert space. At each time step a new 
entangled pair is produced in the above state, and the Hilbert space keeps increasing in size. By 
the end of the evaporation process, on the nth step, we have added 2n qubits to the initial n 
qubits for a total 3n qubits. Since the black hole has completely evaporated and there are only the 
n physical qubits of radiation, the remaining 2n hatted qubits should be interpreted as auxiliary 
degrees of freedom as in the above discussion. Since presumeably the total number of physical 
degrees of freedom should remain fixed at n qubits, at the ith step we have 2i auxiliary (hatted) 
qubits, but the semiclassical analysis does not make any clear identification of the auxiliary qubits 
at intermediate stages in the evaporation. It is clear, however, that by the end of the evaporation 
process all of the black hole (hatted) qubits must be auxiliary Specifying the auxiliary subspace (in 
combination with giving the internal dynamics) corresponds to taking into account back reaction 
on the geometry. 

The Hilbert space is illustrated in Figure [TJ where the intermediate state is shown with some 
highlighted auxiliary degrees of freedom. In the figure, the highlighted region contains some of 
the inital matter qubits (circles) and some of the new inf ailing qubits (squares); this represents 
one possibility. One could also consider cases where for the first n/2 steps only the inital matter 
is auxiliary for example. The parameterization of the auxiliary degrees of freedom at intermedi- 
ate steps should be considered part of a model, so that one may trace out the auxiliary degrees 
of freedom to arive at a fixed-dimensional Hilbert space description. Since we are mostly inter- 
ested in the final state, where the auxiliary space is unambiguous, we do not always specify the 
auxiliary degrees of freedom. In the cases where the evolution is unitary, it should be clear what 
the auxiliary degrees of freedom are. Although nothing profoundly new has been said here, we 
hope this discussion may help clarify potential confusion^] regarding |H0lll5l , and other investi- 
gations IITTtiT3l. 

II.3. The General Model 

We are now ready to present the general class of models that we consider. We start with n 
hatted qubits, and at each step we add a hatted qubit and an unhatted qubit. The number of 
qubits at each step is summarized in Table |Ij 

We model the evolution in two steps: a creation step effected by operators Cf, and an internal 
evolution step effected by Ui acting on the hatted qubits and U{ acting on the unhatted radiation 
qubits. Basis vectors at each step look like 

{\q1q2 ■ ■ ■ q n +i) \qiQi-i ■ ■ ■ qi)} — ^ {\q1q2 ■ ■ ■ q n +iq n +i+i) |gt+igi?t-i • • • qi)} 

Ui ® Uz ) \u\q 1 q2---q n+ iq n+ i +1 }U\qi + iqiq i -i---qi)\ . (2.15) 



6 If we want the initial state to model the matter just before it collapses into a black hole this may not be a good 
assumption, since in a suitably fine-grained description the number of degrees of freedom available to ordinary 
matter is parametrically smaller than the entropy of the black hole |18 19 [. We do not concern ourselves with this 
issue here, but the author is grateful to S. Giddings for pointing this out. 

7 As a specific example, Reference |11| raises the issue of an ever-enlarging Hilbert space as potential issue in the 
analysis of |10|. 



7 



step 


no. BH qubits 


no. rad. qubits 


total no. qubits 


no. 


aux. qubits 


state 





n 





n 







IV>o) 


1 


n+1 


1 


n + 2 




2 


l^i) 


2 


n + 2 


2 


n + 4 




4 


It/>2> 
\ Y*f 


i 


n + i 


i 


n + 2i 




2i 


\i>i) 


n 


2n 


n 


3n 




2n 


\i>n) 



TABLE I. Here we outline the discrete steps in our models. At the Oth or initial step there are n black hole 
(BH) qubits and no radiation qubits. At each step in the evolution, the state is given by the ket in an 
enlarging Hilbert space. 

Of course, one can combine Ci and Ui®Ui into a single operator, but it is useful to break up the 
evolution in this way. Also, there is some physical motivation for thinking about the evolution in 
this way, since the pair creation time scale is roughly ~ M , the black hole mass, while there are 
some conjectures that the internal dynamics of the black hole should be as fast J20j [2T|. (For a 3+1 
dimensional Schwarzschild black hole, the scrambling time is speculated to be ~ M log M, but 
the evolution time step is ~ M.) 

For the majority of the discussion, we focus on the Q and are content to set Ui = Ui = I. What 
properties should the Ci satisfy? We want Ci to preserve the norm and be linear, which means 
that we should require 

(CtfG = I, (2.16) 

where there is no sum on i. Note that this is not the same as Ci{Ci)^ since the Ci have nonsquare 
matrix representations; the above requirement makes Ci an isometric, but nonunitary mapping. 
We also assume that the Ci act only on the hatted black hole qubits and not on the unhatted 
radiation qubits which are far away from the pair creation site. 
We can write the Ci in the following form 

d = \ipx) ® A + \y 2 ) ®P2 + |p 3 ) ® h + |P4> ® A, (2.17) 

where | ipj ) are an orthonormal basis for the created pair qubits, and the P's are linear operators 
which act on the hatted qubits (with implicit i dependence). Following |l2j[l3l, we use the basis 

Wi) = -/=( |6n+i+i) |o i+ i) + |i n +i+i) ) 

l</4) = |0n+i+l) |0 l+ l) - lln+i+l) |lt+l) ) (2.18) 

1^3) = |0n+i+l) 

= |ln+i+l) |0i+l) 1 



for the newly created pair. The constraint in Equation (2.161 implies the following condition on 
fhePs: 

(Ctfd = P}A + P\P 2 + Plh + P\Pa = i- (2.19) 

Note that this defines the Ps as a set of generalized measurement operators acting on the black 
hole Hilbert space. 

A fully specified model, then, entails 
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1. A set of Ps at each step % that satisfy the completeness relation ( 2.19| . 

2. The unitary operators U{ and Ui for each i. 

3. A clear delineation of the auxiliary subspace at each step i. 

The last item is frequently omitted in our discussion; it should be clear for unitary models, and it 
does not make significant differences for the nonunitary models. If one wants to acquire the fixed- 
dimensional Hilbert space description, however, then one must trace out the auxiliary degrees of 
freedom at each step. This gives a very general model space that makes it easy to compare and 
contrast different models of evolution. 



II.4. Physical Motivations 

At this point, the class of models introduced above may seem fairly abstract with little contact 
with the original black hole problem. Let us review the physical motivations for this type of model 
as laid out in [10 j and in ffl2||. 

We consider an initial configuration of spherically symmetric matter that forms a black hole, 
which we expect should be well described by the Schwarzschild solution. The model is based 
on the semiclassical evolution of fields in the background of such a Schwarzschild black hole. In 
order to give a Hilbert space description of evolution, it is necessary to specify a spacelike slicing 
of the geometry so that we can specify the quantum state of fields on each slice. It is important to 
the arguments advanced in [10] that there exists a "nice slicing" of the black hole geometry [22J. 
This slicing avoids the geometry's strong curvature, has sub-Planck scale extrinsic and intrinsic 
curvature, and yet cuts through the initial matter, horizon, and outgoing Hawking radiation in a 
smooth way |Tl0ll22| . Thus, all quantum gravity effects seem to be under control. 

Our model should be considered an effective description of the dynamics on the slicing. As is 
well known, in the presence of curved backgrounds the quantum field theory notion of particle 
becomes observer dependent. If we expand our fields on the slice into modes inside the horizon 
and outside the horizon, then one finds that pairs of particles are created inside and outside of the 
horizon. More explicitly there is a Bogoliubov transformation such that the in vacuum evolves to 
a state of the form exp(7aj nside a^ utside ) acting on the out vacuum. We reduce our problem to an 
essentially two-dimensional one by expanding the modes in spherical harmonics. From a two- 
dimensional perspective, each harmonic corresponds to a different field. Then following ITTI I23II , 
as emphasized in this context in [12], we can use a set of modes that are localized wavepackets so 
that we can talk about locality. This is implicit in the discussion of [ 10]. Moreover, we can truncate 
the Fock space to occupation numbers zero or one. We then are effectively left with a discussion 
of qubits, with |0) representing no excitation and |1) representing an excitation. 

In this description, a pair of particles are created roughly every M in Planck units, with the 
outgoing particles traveling freely outward on the slices and the ingoing particles traveling freely 
inward toward the initial matter that is very far away on each slice. The pair of particles are 
created entangled, as should be clear from the above exponential. This all suggests an effective, 
discrete time evolution with |Tl0llT2"HT5) 

p 1 = I p 2)M = U = U = I. (2.20) 

Because the particles are well-separated on the slice, we do not expect interparticle interactions to 
be significant which is represented by the choice U = U = I. We call this point in model space, 
the Hawking model. Location of the particles on the slice can then be read in the following way 
from the states. Consider, for illustrative purposes, n = 3 with a state of the form 

I9l92g3>initial l 1 405)mfaUing I Oa 1 ! Outgoing ■ ( 2 - 21 ) 
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The first three hatted qubits represent the initial infalling matter. Note that in the Hawking model 
this matter plays no role in the evolution. In general, we imagine that we find some qubit descrip- 
tion of the initial matter, the details of which are irrelevant to our concerns here. The next two 
hatted qubits represent the infalling Hawking radiation as it travels inward on the slices. We see 
that on the first time step, a particle was emitted but not on the second time step. The two unhat- 
ted bits represent the outgoing radiation, so the above implies an outgoing particle was emitted 
on the first time step and then no particle on the second time step. We have written the qubits 
in the above order so that reading from left to right loosely corresponds to traveling outward on 
the slice. This allows us to talk about a coarse form of locality fT0llT2| . Note that in the Hawking 
model the above state would be superposed with several other direct product states. 

As discussed in the Introduction and for this context in |Tl0llT2l|T4l|T5| , * ne arj ove semiclassical 
description of Hawking evaporation is incomplete. In particular, one expects quantum gravity 
effects, backreaction, and interactions to play a role. Because of the nature of the nice slicing 
and the low curvature at the horizon, however, one generally expects all of these corrections to 
be small. That is to say, on the pair creation time scale, one expects from naive estimates that 
the dynamics be e away from the above. This expectation is in considerable tension with the 
expectation that the dynamics are unitary. For instance, we might introduce a set of Ps that act 
on the last emitted ingoing particle. This would suggest that the horizon is not effectively the 
vacuum and still "remembers" the previous emission. One might also allow some mild nonlocal 
interactions inside the black hole via some nearest neighbor UiS. Or, motivated by holography and 
fast scrambling [20ll2Tll24| , one might consider general f/jS, in which case one gives up all notions 
of locality on the slice inside the black hole. The distinction between the initial matter and the 
infalling particles is also lost. Allowi ng g eneral internal dynamics does not affect the argument 



of fiOl and its generalization in Section IV which only relies on the pair creation, P^s, being close to 
the Hawking model. We will always use U = I, since there is no physical motivation to consider 
strong interactions among the outgoing radiation. (Such corrections would also be irrelevant to 
Mathur's bound on the entanglement entropy.) This should become clear after we examine some 
examples. 



III. EXAMPLES 

In this section, we highlight some special points in model space that may be of interest and/ or 
were discussed in the recent literature. One of the main results of this paper is the one-parameter 



family of models presented in Section VI that continuously deforms the Hawking pair production 



in Section III. 1| into the unitary evolution in Section III.3 however, it is also useful to write the 



various models in a common form, so that the similarities and differences are manifest. This is 
especially true when one wants to compare unitary models to nonunitary models. We start with 
the canonical Hawking evaporation model. This should be thought of as the baseline model to 
which all other models should be compared. 



III.l. The Hawking Model 

The standard Hawking evaporation corresponds to creating a new pair in the state \ip\) irre- 



spective of the state of the system, as discussed at length in Section II.4 Thus, it can be written 

as 

Cf = \ip 1 )^I U = U = I, (3.1) 
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and so we can write the Ps as 



Pi=I 



P. 



2,3,4 



(3.2) 



In IHOH , Mathur showed that if the created pair is at most e away from \<px), then the entanglement 
entropy of the radiation continues to grow with each step and thus the final state is mixed and 
unitarity is lost. In this language, the bound shows that if the Ps are small deformations from the 
above, then the final state will be mixed and the evolution will not be unitary. In the sequel, we 
demonstrate this quite explicitly. 

III.2. A Burning Paper Model 

Here, we present a unitary "burning paper" model that is equivalent to the one given in [13jJ^] 

1 



Pl = P2 
P 3 = 



-)=|66) (66| + ^|i6)(i6|-^|io) (oi| 



n+i— l,n+i 



|io) + -)= |66) (io| + -)= |66) (6i| 
Pa = o, 



(3.3) 



n+i—l,n+i 



where the subscript on the brackets indicates that the qubits referred to are the {n + i — l)th and 
the (n + i)th qubits. It should be clear that this model is quite far from the Hawking model. This 
models the creation step as 



Ci = \0 



n+i+l) |li 



|i6)(ii| + -^|66)(io| + -^|66) (6i| 



+ |o. 



n+i+l) |0j 



1 



n+i— 1, n+i 
1 



|()0) (0l)j+^|10)(10|-^|10> (01| 



(3.4) 



n+i—l,n+i 



The important property of the evolution to note is that the (n + i + l)th black hole qubit is always 
0, as is the (n + i)th qubit. Thus these two qubits are "zeroed," and effectively deactivated. We 
use the word zeroed in this sense, even if the qubit under discussion is deactivated to a different 
value. (It could even be something like i mod 2. In this situation, the information is sometimes 
said to be "bleached" out of the state.) It is clear, then, that these two qubits should be thought of 
as the auxiliary qubits at intermediate steps. We'll discuss this a bit more in Section [V] 

We also introduce some interesting internal dynamics. First, we need to move the auxiliary 
qubits out of the way, so that they don't affect the next radiation step. So we first cyclically shift 
all of the qubits 2 positions to the right, thus shoving the two 6s to the two leftmost positions, 1 
and 2. Then, we introduce some dynamics for the physical degrees of freedom. We cyclically shift 
only the nonauxillary qubits to the right by one unit. This defines U so that the model agrees with 
the burning paper model studied in [13J. If we chop off the zeroed qubits, we recover the model 
exactly. 

The first model in |(T2| is in the same class of models. It too zeroes two qubits, one of which is 
the newly created black hole qubit. The main difference being that instead of the radiation being 



Throughout the discussion, the reader may assume that the identity acts on any subspaces which are not explicitly 
shown. 
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determined by the two rightmost hatted qubits, it is instead determined by the leftmost qubit. 
This model may be written as 

A=A = -s|6i+i)(6 i+ i|®ti 

V2 (3.5) 

P3 = |6j+i) (lj+i| <8> u', 

where u and u' are unitary operators acting on the remaining hatted qubits. While it is not stated 
in fl2fl, we should require that u and u do not mix the first i + 1 or the last i auxiliary hatted 
qubits with the remaining physical qubits. The simplest case is to take u = v! = I. In this model, 
the entanglement entropy of the radiation is always zero, which contrasts with our expectations 
from 0. 



III.3. "Nonlocal" Unitary Evolution 

In fl2j , Giddings presents three unitary models of evolution. We focus on the second model 
that he presents. This second model can be written in our notation as 



-Pi — |02i+l02i+2) (02i+l02i+2| 

Pi = |66> (ii| 
P 3 = 1 66) (6i| 
A = |66) (io| 



U = I, (3.6) 



where we have suppressed the (2i + 1) and (2i + 2) subscripts in all but the first P. One sees that 



as in the models presented in Section III.2 two hatted qubits are zeroed at each step. In this case, 
they are the (2i + l)th and (2i + 2)th qubits. Thus as the evolution progresses the hatted qubits 
are gradually put into a fiducial form. By the ith step, the first 2i qubits are zeroed, and should 
be thought of as auxiliary Note that the above evolution rule breaks down on the penultimate 
step, when there is only one nonauxillary qubit left. By then, we expect the black hole to be on the 
Planck scale, and so we can just emit the last qubit freely. This model corresponds to 9 = ^ in the 



model presented in Section VI 



This model is nonlocal when one considers the model in the original nice slicing of the black 
hole. In this context, the (2i + l)th and (2i + 2)th qubit are very far from the pair creation site 
at the horizon. Note that this property is shared by the model in Equation ( [3.5^ and to a lesser 
extent the model in Equation The difference being how far away the zeroed qubits are. 

One can either interpret these unitary models as nonlocal interactions transmitting information 
far down the nice slice to the horizon [12J, or in terms of fuzzball microstates altering the state at 
the horizon, or as burning paper; these information theoretic models are too crude to distinguish. 
This point is discussed futher in the Conclusion. 

III.4. A Pure, but Not Invertible Model 

There are several of models of this kind. The simplest to consider is 

A, 2 ,4 = o P 3 = i u = i. (3.7) 



In this model, regardless of the state of the system, the new pair is created in the state |<^3). The 
state 1993) is not entangled, and thus we can think of this model as zeroing the new black hole 
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qubit and the new radiation qubit. Instead of putting the internal qubits into a fiducial form, we 
put the radiation into a fiducial form. We hope that this convincingly demonstrates that purity of 
the final state does not ensure unitarity. 

Another interesting example of (almost) pure but not invertible evolution is the model in Equa- 
tion \3.5) , with 

M+2,n+i /o Q s 

u = u = b 1 , (3.0) 

where Sl +2 ' n+1 is the operator that cyclically shifts the (i + 2)th through (n + i)th qubits to the left. 
For example, consider the evolution: 

|6<?1<?2 • • • Qn-20) 

^|60gig 2 ---<7 n _ 2 6>|0) 

% 1000^2 • ■ • <? n _ 2 0) |00) (3.9) 
^> 10000^1^2 • - - g„_ 2 6> |000) 



In the above example, the radiation is never entangled with the black hole degrees of freedom, 
but its state is only determined by the first and last qubit of the initial state. Thus, the evolution 
is not unitary. This is a potential problem with the model even as it is defined in fl2"|. What 
happened? In essence, the model keeps "reading" and zeroing the same qubit, which has the net 
effect of zeroing the radiation qubits as well as the new black hole qubits. This illustrates that 
unitary evolution can be lost when auxiliary qubits mix with nonauxillary qubits. 
This model is not quite pure for arbitrary initial states, since 



Co 
Ci_ 

C 2 . 



|l<7l<72 • • -<?n-2l) 

\6iq 1 q 2 ---q n -20) |1> 

|66o<M 2 ---(z n _ 2 6)|ii) , (3.io) 
|6666<M2 • • • q n -20) |oii) 



which gives a pure final state, but if one considers a nontrivial superposition of the above initial 
state and the one in ( |3.9| then the final state is mixed. Note that the radiation only carries two 
qubits of information about the initial state, and so the entropy of the final state is very small 
although nonvanishing on some initial states. It is in this sense that we call the evolution almost 
pure. 



III.5. An Impure Model 



A generic model that one constructs leads to impure evolution. When thin king about the 



different forms that Ps can take, a particularly natural variation on the Section 
consider might be 

A = |66> <66| 
A = 

A = |6i) (6i| 
A = |i6)(i6|, 



III.3 



model to 



(3.11) 
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where as in (3^1 the above operators act on the (2i + l)th and (2i + 2)th qubits. It should be clear 
that this model leads to mixed states. Note that it is also a large deformation from the Hawking 
model. It is easy to see that this model is not invertible either, by considering 1 00) and |11) as 
initial states. 



III.6. Mathur-Plumberg Shift-Anti-Shift Models 

In [[13"), Mathur and Plumberg present several models. We can write their "Model A," in the 
following way. Let Tj be the operator that cyclically shifts only the newly created (not the first n) 
hatted qubits to the right by j. Then, the model is of the form 

A = Aiii A = A 2 T_i A,4 = 0, (3.12) 

where the completeness relation requires 

|Ai| 2 + |A 2 | 2 = l. (3.13) 

More generally, we can replace T\ and T_i by other unitary transformations that act on the hatted 
qubits. The smallness of the corrections to Hawking is loosely determined by A 2 ; see [13j for more 
details and numerical results. 

The "Model B" of |l3l may be written as 

A = Aiii® |i i+ i)(i m | + J<g>|6 m )(6 i+ i| A = A 2 f_i<g)|i i+ i)(i m | A, 4 = o, (3.14) 

where one can confirm that the completeness relation imposes the same constraint ( 3.13} . 

Note that neither of the above models zeroes any qubits, and so one does not expect informa- 
tion about the initial matter to be transmitted out in the radiation. 



III.7. Mathur "Ising" Model 



In |Tl5l , Mathur presents a model that can be mapped onto the one-dimensional Ising model 
and thus solved analytically. The model can be written in the form 



A = Ax [|o> (6| + |i> {l\] n+i A = A 2 [|6> <6| - |i) {i\] n+t 

where Ai and A2 are related to the parameters a and b in jT5] via 

g a _j_ ^ Q sfi 



3,4 



0, 



(3.15) 



Ai 



A, 



The above rule is not valid for the first step, for which we use the Hawking model (Ai = 1 and 
A 2 = 0). The parameters Ai and A 2 are subject to the same constraint as before 



|Ai| 2 + |A 2 | 2 = l. 
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(3.16) 

both set Pi and P 2 to 



In fact, this model is qualitatively the same as the model in Section 
unitary transformations, and P3 = P4 = 0. The expression for the final state entropy can be 
written in the form [15 j 



S = nlog2 - (n- 1) 



ae 2a + be 2b 



(3.17) 



14 



Let us note that for no value of a and b that solves the completeness relation does this vanish. 
Thus, no matter how large the corrections are in this model, unitarity is lost; however, one can set 
Ai = A2 in which case the entanglement entropy of the radiation is log 2 for all time. In this case, 
the pair creation (after the first step) is given by 

Ci = |6 n+ j + i) |0i) tg> |6 n+i ) (6 n+ j| + |i n +i+i) |ii) ® |in+i) (in+i| • (3.18) 

One sees that after the first step, no further entanglement between the hatted and unhatted qubits 
is generated. One might think there would be more entanglement generated when the above acts 
on qubits which are in a superposition of 1 and 0, but since the above only acts on previously 
created pairs we do not have to worry about that issue. One can imagine that the model breaks 
down when the black hole becomes very small and the last qubit is emitted freely, so this model 
is effectively pure in this case. Let us note that keeping the entanglement entropy of the radiation 
constant is in stark contrast with expectations we have from Page [8J; there is no characteristic 
rise and fall of the entanglement entropy of the radiation. Furthermore, the final state consisting 
entirely of radiation is completely independent of the initial matter that formed the black hole; the 



evolution is (almost) pure but far from invertible like the model discussed in Section III.4 This 
model, while interesting to consider, never leads to unitary evolution, even when one considers 
arbitrarily large deformations from the Hawking point in model space. 



IV. REVIEW AND GENERALIZATION OF MATHUR'S ARGUMENT 

In Pol , Mathur argues that small corrections to the pair creation process are insufficient to 
restore unitarity. More specifically, he demonstrates that small corrections don't accumulate. Let 
us define S{ as the entanglement entropy of the radiation with the rest of the system at step i. 
If the black hole evaporation process is unitary, then in the limit of large n we expect Si to rise 
linearly with i until about the halfway point, i = n/2, and then rapidly turn over and fall linearly 
to zero on the final step [9j. On the other hand, if one uses the Hawking model of evolution C H , 
then one sees that Si increases by log 2 at each time step: 

^Hawking = . lQg2; (41) 

there is no turnover. This is also the maximum entropy that is possible for the i radiation qubits, 
which indicates that the radiation carries no information about the initial state. 

The central insight in Mathur 's argument ffT0| is that the marginal increase in entanglement 
entropy, 

ASi = S l+1 - S u (4.2) 

varies smoothly with small deformations away from the Hawking model, and thus a large defor- 
mation is needed to make AS negative if one starts from the Hawking model's log 2. In IfTUl , only a 
-fVtype deformation to the Hawking model was explicitly considered; however, we demonstrate 
below that the the argument generalizes to all deformations in the model space. 
In particular, we claim in the class of models discussed in this paper, if 

\\d -Cf || < e < 1, (4.3) 

then 



ASi > log 2- k £ , 



(4.4) 
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where k £ is parametrically small, positive, and vanishes as e goes to zero. In fact, we demonstrate 
below that k £ behaves as no worse than k £ ~ — 9eloge as e approaches zero. Above, || • || is the 
operator norm. For an operator O, \\0\\ is the square root of the largest eigenvalue of O^Oj^] 

The proof follows that in [ 10 1 with some minor modifications. To begin, we use strong subad- 
ditivity of (von Neumann) entanglement entropy. Throughout, without loss of generality, we set 
U = U = I since they do not affect the entanglement entropies. Let us consider the (i + l)th state 

= Ci \tl>i) (4.5) 

Let Ri denote the first i emitted radiation qubits, r denote the (i + l)th emitted radiation qubit, 
Bi denote the first n + % black hole qubits, and b denote the (n + i + l)th black hole qubit. Then, 
strong subadditivity implies 

S(R l U r) + Sir U b) > S{Ri) + S{b) . (4.6) 

By definition, S(Ri U r) is simply Si+i, whereas since the Ci acts trivially on emitted radiation 
S(Ri) = S{. Thus, we can write the above as 

ASi> S(b) - S{rUb). (4.7) 

Note that for the Hawking model S(b) = log 2 and S(r U b) = 0, and the bound is saturated. 

Now, we need to use the condition on Ci to place bounds on S(b) and S(r U b). The operator 
norm is compatible with the Hilbert space norm, which means 

\\(C l -C H )\X)\\ < ||a-C|||||A>|| (4.8) 

for all | A). Applying this to \ipi) gives the condition 

|||^ +1 )-Cf |^)|| <s. (4.9) 

We may write the two kets in the form 

= at |Ai) + q 2 \(p2) |A 2 > + as |^3> |A 3 ) + a 4 |^4> |A 4 > , Cf \^) = \(p x ) |A ) , (4.10) 

where the | A, ) are normalized, but not necessarily orthogonal kets in the n+2i qubit space R4 U Bi . 
Of course normalization demands that 

Kl 2 + M 2 + 1«3| 2 + 1«4| 2 = i- (4.11) 

One can show that the condition ( |4.9| implies that 

Re(ai(Ao|A a >) > l-j, (4.12) 

and since | (Ao| Ai) | < 1, we see 

1 - y < \ai\ < 1. (4.13) 

This, in turn, implies that 



1 - 5 2 < \a\\ z < 1, \a 2 \, \a 3 \, \a^\ < S = e-i/1 - (4.14) 



e 2 



9 In fact, the result is unchanged by using any other norm that is compatible with the Hilbert space norm. 
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where we have defined 5 for convenience. Note that 5 is what should be directly compared with 

e in on. 

Let us use the above to place an upper bound on S(r U b). The corresponding reduced density 
matrix is given by 



Prb 



\ai\ 



a\a\ (A 2 |Ai 



«2 



a x a\ (A 3 |Ai) 
a 2 a* 3 (A 3 |A 2 ) 



a^a 2 (Ai|A 2 ) 

a*a 3 (Ai|A 3 ) a*a 3 (A 2 |A 3 
\a^a4(Ai|A4) aJja^ (A2IA4) aga4(A 3 |A4) 



\ a 3\ 



aial (A 4 |Ai) N \ 
a 2 a| (A 4 |A 2 
a 3 a| (A 4 |A 3 



a 4 



(4.15) 



/ 



We can now use Audenaert's optimal generalization H25l of Fannes' inequality which places a 
bound on the difference in entropy of two <i-dimensional density matrices p and uj^Let T be the 
trace distance between p and a, then [25 1 



\S(p)-S(a)\<Tlog(d-l 

where the trace distance is defined as 

_ 1 



TlogT-(l-T)log(l-r), 



tr 



^{p-cr)Hp 



(4.16) 



(4.17) 



or one-half the sum of the absolute value of the eigenvalues of p—a. Note that the above definition 
of the trace distance differs by a factor of 2 from some references. With the above normalization 
< T < 1 for all unit-trace density matrices p and a. In this case we consider the a to be the 
density matrix with |aj.| = 1, and p to be p r b- Next, we may use Gershgorin's circle theorem to 
bound the eigenvalues, A, and therefore the trace distance, T. Gershgorin's theorem tells us all the 
eigenvalues must lie in the union of discs in the complex plane centered on the diagonal entries 
with radii given by the sum of the absolute value of off-diagonal entries for each row. For instance 
the first row gives a disc 



Di : |A - (|ai| 2 - 1)| < |ai| \a 2 \ \ (A 2 |Ai) | + | ai| |a 3 | | (A 3 |Ai) | + |ai| \a^\ \ (A4IA1) 



(4.18) 



Applying our inequalities on the components and the hermiticity of p — a (and therefore reality 
of its spectrum), we can find an interval that must include any eigenvalues in the above disc: 

I 1 = (-35 - 5 2 ,35). (4.19) 

One finds the remaining three rows can be encompassed by the interval 

I 2 = ( y S- 25 2 ,5 + 35 2 ), (4.20) 

and so all eigenvalues must satisfy 

\X\< 38 + 5 2 . (4.21) 
Thus, we may conclude that the trace distance must satisfy 

T < 2(36 + 5 2 ), (4.22) 

and therefore 

S(r U b) < -xi log (ixi) - (1 - xi) log(l - x\) x x = min (§, 2(35 + <5 2 )). (4.23) 



One could instead use Fannes' inequality for a weaker bound with stronger restrictions on e. 



17 



The 3/4 comes by finding critical point of the right-hand side of (4.16 1 as a function of T. 

The bound on S(b) can be derived in analogous fashion. The state may be written out 

as 

IVi+i) = |6)|xo) + |i) 

The reduced density matrix can be written as 



Ixo; 



Xl 



Ixi) 

|0)|A 2 ) + a 3 |l) 



I A, 



(4.24) 



|l)|A 2 ) + a 4 |0)|A 4 ) 



Pb 



(Xolxo) 
(Xolxi) 



(Xilxo) 

(xilxi) 



(4.25) 



As before we can place bounds on the above components 



(Xolxo) 



(xilxi) 



< 5 + 



<5 + 



2 

S 2 



(4.26) 



( X i\xo)\<V2(5 + 6 2 



One can once again use the Fannes-Audenaert inequality along with Gershgorin's theorem to 
bound S(b), where a is the a\ = 1 density matrix, 1/2. One finds the trace distance satisfies 

T< (l + y/2)6 + {\ + V2)5 2 , (4.27) 

and this gives 

S{b) > log 2 + x 2 logx 2 + (1 - x 2 ) log(l - x 2 ) x 2 = min (§, (1 + y/2)6 + (± + V2)S 2 ) . (4.28) 

Finally this allows us to write 

k £ = -xi \og(\xi) - x 2 log x 2 - (1 - xi) log(l - xi) - (1 - x 2 ) log(l - x 2 ), (4.29) 

where recall x\ is defined in Equation (4.231, x 2 in Equation (428), and 5 in Equation (4.14|. For 
asymptotically small e, k £ ~ — (7 + v / 2)eloge; and in fact for small but finite e, k £ < — 9eloge. 
Numerically, one finds that k £ first surpasses log 2, thus allowing the entanglement entropy to de- 
crease for e .02. Furthermore, one finds k £ reaches 2 log 2, thus allowing the maximal marginal 
decrease of entanglement entropy for e .05. For larger e, the inequality with k £ given in ( 4.29| 
places no restriction on the marginal change in entanglement. 

Since we are not especially interested in making the tightest possible bound or even the above 
numerical values, we may as well write the bound in slightly less unwieldy form, 

ASi > log 2 + 9e log e e < 1 . (4.30) 

This establishes the claim. Let us note that the above bound's asymptotic behavior is weaker than 
the inequality derived in Jlpj^jas a consequence of using more general arguments to include 
arbitrary perturbations. On the other hand, the result presented here is stronger in the sense that 
iTlOl finds only a leading order result valid to order 0(e 2 ) whereas the bound ( |4.4| | with k £ given 
in (429) is valid for finite e € (0, 1). The above bounds could possibly be strengthened with more 
workj 12 |however, that is irrelevant to the basic claim that small corrections to the low energy pair 
creation process cannot restore unitarity. 

11 An equivalently strong bound here would be k E = 28. 

12 For example, one might make progress by direct computation of S(r U b) and S(b) as was performed in ] 10 \, but this 
would involve solving an eigenvalue problem for a four-dimensional density matrix with arbitrary coefficients. 
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V. REQUIREMENTS FOR UNITARITY 

While it is interesting to think about the different kinds of evolution that one could have, per- 
haps the most interesting question to ask is what kinds of models are unitary or, equivalently 
what sorts of corrections to the Hawking evolution can restore unitarity. Above we see that small 
corrections to the evolution cannot restore unitarity; this gives a necessary condition that the cor- 
rections are large. It would be nice to also have some sufficient conditions, since it is clear that not 
every large correction one could consider leads to unitary evolution. 

In order for the evolution to be pure, we need the final state (including the auxiliary qubits) to 
be a direct product of the form 

IV>n> = \4>) ® |X) , (5.1) 

where \cf>) is a state in the 2n-qubit auxiliary space and |x) is the state of the physical radiation 
qubits. Our first observation is that \4>) should be independent of the initial state. Suppose that 
this were not true: 

| W) |0(D) s | x (i) } 
\4&) » \4> (2) ) ® lx (2) >' 

which seems fine until one considers an initial state which is a superposition of the above two 
states; the final state is then mixed when one traces out the hatted qubits. (This argument assumes 
that the \x)s are linearly independent so that the evolution is invertible.) This is basically a variant 
of the no-cloning theorem. 

The next observation is that we need the final radiation state | \) to be a unitary transformation 
of the initial state \ifio)- Let F be the total map from initial state to the final state, then F is a linear, 
isometric (norm-preserving) mapping from n hatted qubits to 2n hatted plus n unhatted qubits. 
From the above, unitarity demands that 

^unitary : |$ -> \4>) ® \x) \x) = U \$) , (5.3) 

where U in the above is a unitary transformation from the initial n hatted qubits to the final n 
unhatted qubits, and \(f>) is fixed. The total map F is just the product of all the CjS and U ® Us. 
All of the hatted qubits have to be zeroed or bleached, and the information stored in the initial 
matter transferred to the radiation. 

Let us think about how we can zero or bleach the hatted qubits. We have n steps to project 2n 
qubits to a unique state with the P l s; the unitary Us clearly cannot zero qubits. At each step, we 
can zero at most two qubits. If C bleaches some subspace to a state \a), then it may be written as 

C=\a)®0, (5.4) 

for an unspecified operator O. If \a) is a p-qubit subspace, then O maps n + i qubits ton + i + 2—p 
qubits and must satisfy 

] = /; (5.5) 

this is only possible ifn + i + 2— p > n + i, immediately implying p < 2. 

Our key observation is that the desire to zero the hatted qubits is in tension with the complete- 



ness relation (2.19 (. Since we only have four Ps, at any given step the best we can do is project out 
a four-dimensional subspace, or two qubits. It is this tension that connects the need to zero qubits 
with the requirement to have large corrections to the Hawking model. This might help elucidate 
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the results in [10 j. Note that this is very much in agreement with the picture presented in FigurejTJ 
wherein at each stage there are two new auxiliary qubits. For the state to be pure, these auxiliary 
qubits must be zeroed. 

Let us note that there are two different ways to zero two qubits at each step, although the dis- 
tinction is not a ctually that sharp when one considers the full CjS. In the burning paper model 
we use only three Ps, which zero one qubit. The three Ps were chosen to ensure 
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of Section 

that the newly created q is also zeroed. The second way is illustrated in Section III.3 in which all 
four Ps are used to zero two old qs. One can consider various unitary transformations, however, 
these are the only two qualitative kinds of models that lead to a pure radiation final state. Re- 
member that which Ps get used tell us which pair state is created at the horizon. It is impossible 
to preferentially use only P\ and simultaneously have unitary evolution. 

As we saw in Section III.4 it is possible for the evolution to be pure, but not invertible. In 
the model in Equation ( |3.8| l, qubits that were zeroed in previous steps mixed with nonzeroed 
qubits. When we zero the qubits, we are then thinking of them as auxiliary degrees of freedom 
that should be erased in the operator-sum description ( |2.6| |. Thus, it does not make physical sense 
to allow mixing with the auxiliary degrees of freedom if one wants unitary evolution. 

The requirements outlined above for purity and invertibilty should ensure unitary evolution. 



VI. A ONE-PARAMETER INTERPOLATING MODEL 



We can interpolate between the Hawking model in Section III.l and the unitary model in Sec- 
tion HlO] via 



Pi = cos 91 + (1 - cos I 
P 2 = sin<9|66) 

p 3 = sin e 1 66) (6i| 

P 4 = sin0|66> (io| 
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2i+l02i+2) (02i+l02i+2| 



u = i, 



(6.1) 



with 9 = giving Hawking's evolution and 9 = | giving the unitary evolution in Equation (3.6 1. 
(Once again we suppressed subscripts on the qubits after the first line.) We may write 



(C - C H )\C - C H ) = 21 - A - PI = 2(1 - cos 9) (I - |66) (66| ) , 



(6.2) 



so that one finds 



IC-C 



H I 



2 1 sin 



(6.3) 



We clearly see that this is in accord with Mathur's argument and its generalization in Section IV 



One of the main results of this paper is the above model, which continuously connects the 
Hawking model to a unitary model, clearly illustrating that they are far apart in model space. 
Previous efforts to illuminate Mathur's bound IIT3l [T5| , considered different types of small cor- 
rections and showed that they did not significantly affect the entropy of the final state; however, 
they did not consider corrections that when made large would give a unitary "burning paper" 
type model. The above model fills this gap. In Figure |2j we plot the second Renyi entanglement 
entropies of the radiation as a function of 9 for three different initial states. Recall that the second 
Renyi entropy is defined as 



Sa(Awd) = -!ogtr(Pr ed ), 



(6.4) 
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and is a positive-definite measure of entanglement that vanishes if and only if p re .^ is pure. For 
computational purposes, however, it is a bit more convenient than the traditional von Neumann 
entropy. It should also be noted that S2 gives a lower bound for the von Neumann entropy. 
Moreover, when one examines the von Neumann entropy it behaves qualitatively similarly. 




(a) |66006) 



(b) |666ii) 



(c) 



FIG. 2. Here we present the second Renyi entropies of the unhatted radiation qubits as a function of 8 for 
the model in Equation 16.1 1 with three different initial states. Intermediate steps are dashed curves and the 
final step is solid; the steps are shown in the color order (red, yellow, green, blue, purple, black). The point 
6 = corresponds to the canonical Hawking evolution, while 6 = § corresponds to the unitary model in 
Equation ||3.6). 



In Figure|2j one sees that at 9 = (Hawking model) the entropy rises by one at each step, except 
for the final step where it drops down by one. This is an artifact of the way we end the evolution, 
since the model breaks down on the penultimate step. We chose to just emit the last qubit freely, 
which makes sense on physical grounds since the black hole should be quite small by that point; 
one can easily imagine large corrections in the final stage(s) of black hole evaporation. At 9 = 
\, we have the unitary model, and the entanglement entropy of the radiation has the expected 
rise and fall. For no 9 close to 9 = 0, however, does the entanglement entropy of the radiation 
fail to increase (excluding the final step). Of course, a 5 qubit initial state is not realistic for the 
macroscopically large black holes we are thinking about; however, it suffices to demonstrate the 
qualitative behaviour which should not change as n is increased. The author was limited by the 
computational power required, which grows quite rapidly with n. 



VII. CONCLUSION 



We have presented a very general framework that provides a natural, unifying language to 
compare information-theoretic models of black hole evaporation. The framework involves de- 
scribing the dynamics of pure state vectors \tp) in an ever enlarging Hilbert space. The only 
constraint on the models is the completeness relation ( 2.19) . In Section II.2 we explain how to 



interpret this model in terms of potentially mixed evolution in a fixed dimensional Hilbert space: 
one must trace out some auxiliary degrees of freedom to arrive at a dynamical equation like in 



Equation (2.6 1. Part of a full model, then, involves specifying the auxiliary degrees of freedom. 
While at intermediate steps this can be ambiguous, by the end of the evolution we are left only 
with radiation and thus any nonradiation degrees of freedom are by default auxiliary This ex- 
cludes remnants or other scenarios where one can identify physical degrees of freedom that the 
Hawking radiation is entangled with at the end of the evaporation. 



In Section III we show how to write a number of interesting models in our unifying nota- 
tion. Many of them had been introduced and studied previously in the literature. These models 
illustrate some of the key obstructions and requirements to have unitary evolution. 
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In Section |Vj we discuss the requirements for unitarity, and summarize a set of sufficient con- 
ditions. A key backdrop to our discussion, and indeed the whole paper, is a recent theorem IflOl 
showing that small corrections to the Hawking model cannot give unitary evolution. A corollary 
is that large corrections are necessary to have unitary evolution. As we show, this is, however, 
not a sufficient condition; unsurprisingly, there are many large corrections that fail to give uni- 
tary evolution. One interesting observation is how the unitary requirement that internal qubits be 
zeroed becomes connected to corrections to the pair creation via the completeness relation ( |2.19[ ). 
This may help elucidate the results in IITOl . 

Finally, in Section VI we give a one-parameter family of models that continuously interpolates 
between the Hawking model and a unitary model. In terms of this parameter, one can clearly see 
that the unitary model is far from the Hawking model, thus illustrating the theorem in IITOll . 

One of the key points emphasized in flOllMl is that the nice slicing of the Schwarzschild solu- 
tion implies at best small corrections to the Hawking model in Equation ( p.lf , and therefore a loss 
of unitarity evolution. To restore unitarity, large corrections are required of the form discussed 
here; however, one must show why these corrections arise in the black hole and not in all of our 
earth-based experiments and observations. It seems quite difficult to do this, since in the nice 
slice construction no geometric quantity is large. The only quantity that seems to be large is the 
number of degrees of freedom or number of particles required to form the black hole, but this is 
not a basic geometric quantity. 

Let us further note that for our discussion there is a factorization of the internal dynamics and 
the pair creation dynamics. The internal black hole dynamics can be as nonlocal, or scramble 
as rapidly as one wants, but whether the evaporation process is unitary or not is (modulo a few 
caveats mentioned in Sections III and|V]l entirely determined by the pair creation process. Thus, 
the discussion in [20} [2111241 [26 1 is not directly relevant to our concerns here, although it is impor- 
tant for better understanding black holes. The pair creation process is localized near the horizon, 
where for a large black hole the geometry suggests one can trust the semiclassical approximation 
even if one might doubt its validity deep within the black hole. On the pair creation time scale, 
however, it is precisely this physics that needs an order unity correction [10j. 

One other issue that one may wish to raise in our discussion is the issue of conservation 
laws ffTTll27l . While we hope we have sufficiently addressed the issue of an expanding Hilbert 
space as raised in [[rT| , one may still be concerned that we haven't discussed conservation of en- 
ergy (or angular momentum, electric charge, etc.). These issues are rebutted in [14J. Let us note, 
however, by not discussing the original spacetime physics and the resulting at most small correc- 
tions, the fundamental issue has been totally elided. While the results presented in iTlTl |27|| are 
interesting unto themselves, they do not provide a plausible physical mechanism to modify the 
pair creation process to get the dynamics they suggest. 

If, as suggested by string theory, or from other considerations, we think black hole evaporation 
is a unitary process; then, the pair creation process must not strictly adhere to the causal structure 
on the Schwarzschild nice slicing. There are two obvious frameworks (ignoring the possibility of 
remnants) to discuss this deviation: fuzzballs or nonlocality. 

The fuzzball proposal (see 12814321 for reviews) suggests that the black hole metric is only an 
effective geometry that approximates e 5BH microstates. The microstates differ from each other on 
the horizon scale, thus large corrections to the Hawking evolution are anticipated and informa- 
tion is transmitted from local excitations. How the fuzzball proposal relates to these qubit models 
is discussed in ITTOI [1314T5I . The main point being that since the geometry at the would-be hori- 
zon depends on the internal state, their is a physical mechanism to get large corrections to the 
pair creation process. Since the fuzzball's interior geometry (and the whole causal structure) is 
quite different from the original black hole solution in which the nice slices were constructed, 
one should probably not interpret them with the original notion of locality for the internal qubits 
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discussed in Section II.4 Moreover, let us note that by adding nontrivial dynamics of the inter- 
nal fuzzball structure via a U to the model in \3.6), one can effectively change which qubits get 



emitted. This is the point referred to at the end of Section III.3 



There is one explicit family of (non-extremal) fuzzball microstates 11331 for which one can un- 
derstand the bulk Hawking emission process. As first suggested in [34], the geometry's ergore- 
gion instability |35| can be interpreted as a Bose enhanced version of the Hawking instability for 
the corresponding black hole. This explanation was justified by comparing gravitational emission 
to the dual CFT emission process for both ergoregion emission from the fuzzballs and Hawking 
radiation from the corresponding black hole fl34l [36 -39] . In [34], a toy model was presented for 
the ergoregion emission, based on the CFT description. The toy model consists of a set of two- 
level atoms that can spontaneously emit or absorb photons. In the geometric description, the 
de-exciting atoms correspond to accumulating particles in the ergoregion that decrease the ge- 
ometry's mass and angular momentum. Loosely, in our language, the toy model is in the class 
discussed in Section III.2 To properly capture the Bose enhancement, however, is a bit trickier. 

In fl"2| , several unitary models of evolution (some of which were discussed here) were pre- 
sented, motivated by proposed nonlocal physics on the Schwarzschild background. As men- 
tioned in BUI , it is unclear what sets the scale of the proposed nonlocality so as to ensure it oper- 
ates in the black hole background but not in everyday low-energy experiments. While the sorts 
of models discussed here remain too crude to distinguish between nonlocal physics or fuzzball 
microstates, their utility lies in their generality, which serves to sharpen our information theoretic 
understanding of black hole evaporation. One obvious task that remains is to translate Mathur's 
bound and its generalization in Section IV into a sharper, quantitative statement about the break- 
down of the semiclassical limit of quantum gravity. Since the entire discussion has been in a 
Hamiltonian framework, it would be especially nice to have analogous bounds on the path inte- 
gral. 
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